Ultradiscrete Soliton Systems and 
Combinatorial Representation Theory 



Reiho Sakamoto 

Department of Physics, Tokyo University of Science, 
Kagurazaka Shinjukuku Tokyo Japan 

1 Introduction 

This lecture note is intended to be a brief introduction to a recent development on the 
interplay between the ultradiscrete (or tropical) soliton systems and the combinatorial 
representation theory. We will concentrate on the simplest cases which admit elementary 
explanations without losing essential ideas of the theory. In particular we give definitions 
for the main constructions corresponding to the vector representation of type A^^K 

This note is organized as follows. In Section [2] we give a definition of the simplest ex- 
ample of the box-ball systems. In Section [3] we explain a relationship between the box-ball 
systems and the crystal bases of the quantum affine algebras. In Section H] we give the 
definition of the rigged configuration bijection for the vector representation of type A^^\ 
In Section E] we see that the rigged configurations give the complete set of the action and 
angle variables for the box-ball systems. This is the fundamental observation in the recent 
development on a relationship between the box-ball systems and the combinatorial repre- 
sentation theory. In Section |6] we explain basic properties of the box-basket-ball systems 
which are recently found generalizations of the box-ball systems. The characteristic prop- 
erty of the system is that it is a mixture of fermions and bosons with mutual interactions. 
Finally, in Section [3 we give comments on generalizations and further developments of 
the materials discussed in this note. 

This is the lecture note prepared for the conference "Algebraic combinatorics related 
to Young diagrams and statistical physics" held at International Institute for Advanced 
Studies (Kyoto) during August 6-10, 2012. The author is grateful to Professor Masao 
Ishikawa, Professor Soichi Okada and Professor Hiroyuki Tagawa for the kind invitation 
to the conference and the warm hospitality. 



2 The box-ball systems 

In this section, let us define the simplest case of the box-ball systems introduced by 
Takahashi and Satsuma |TS] . The box-ball systems are prototypical examples of the 
ultradiscrete soliton systems. Originally the ultradiscrete soliton system is a class of 
discrete dynamical systems obtained by the ultradiscrete (or tropical) limit of the ordinary 
soliton systems |TTMS] . In this article we are interested in ultradiscrete soliton systems 
which admit combinatorial interpretations. 
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Following the box and ball interpretation of the system |T3] , we prepare boxes which 
can accommodate at most one ball within each box. We put many such boxes on a line 
and put finitely many balls of the same kind to the boxes. We regard this configuration 
as the initial state of the system. Then we perform the time evolution of the state by the 
following algorithm. 

The time evolution of the box-ball system: Consider each ball from left 
to right and move the ball to the next available empty box. Each ball is moved 
exactly once. 

Here if necessary we put enough many empty boxes on the right of the given state in order 
to keep the balls within the state. We give an example of such time evolution starting 
from the top row and proceeding downwards. 
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In the box-ball system, we regard consecutive balls as solitary waves. For example, the ini- 
tial state in the above example contains two solitary waves of length 3 and 1, respectively. 
Then our interpretation of the above example is as follows. If there is no interaction 
between waves, they move at velocity equal to each length. In the course of the time 
evolution, solitary waves make collision with each other, though they retain their original 
shapes after the collision except for the changes in the positions compared with the pos- 
sible positions we would have if there is no interaction between waves. Such properties of 
the waves of the box-ball systems are characteristic of the soliton systems (see, e.g., |T4j ) 
and we will call such waves solitons. 

Here we give a short list of remarks on the early papers. During 1980's, there were 
several attempts of finding cellular automata with solitonic properties. A typical example 
of such researches is the filter automata introduced by Park, Steiglitz and Thurston |PST] . 
In 1990, Takahashi and Satsuma |TSj introduced the simplest case of the box-ball systems 
and Takahashi |T3j described the algorithm in terms of the boxes and balls. The above 
definition corresponds to the original Takahashi-Satsuma box ball system. A relationship 
between the Takahashi-Satsuma box-ball system and the ordinary soliton systems includ- 
ing the KdV equation is discovered by Tokihiro, Takahashi, Matsukidaira and Satsuma 
[TTMSj via the limiting procedure called the ultradiscrete limit. A connection with the 
Toda equation is discussed in [NTT]. Such connections between the box-ball system and 
the usual soliton systems show the classical integrability of the box-ball system. 

Takahashi's box and ball algorithm provides several generalizations of the original box- 
ball system. For example, in |T3] an internal degree of freedom for the balls (balls with 
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different colors) is introduced. A connection with the Toda equation in the generalized 
context is discussed in |TNS] . The other degrees of freedom called a carrier |TM] or a 
capacity of boxes |TTM] are also introduced. Such combinatorial interpretations of the 
time evolutions give nice intuition about the models in many cases. 



3 A connection with the crystal bases 

A very important fact IHHIKTTl IFOYj about the box-ball systems is that their dynamics 
is in fact governed by the Kashiwara's crystal bases |K2j for the quantum affine algebras. 
This formalism includes all extensions of the box-ball system which are mentioned in the 
last section. Although the formulation does not depend on the types of the algebra, we 
will concentrate on the simplest case A^^^ here. 

In order to describe the formulation, we need to consider more general boxes which 
have capacities more than one. Let (a, b) represents the box of capacity a + b containing 
b balls. Then the state (a, b) can accommodate extra a balls. Let us denote the set of all 
such states as 

B^'' := {{a, b)\ a, be Z>o, a + b = s} (1) 

which we call crystalsj^ In particular we call B^'^ the crystals for the vector repre- 
sentation. In this coordinate, the states of the box-ball system in the previous section 
are sequences of balls (0,1) and empty places (1,0). Then we represent the states as 
(1, 0) (8) (0, 1) • • • where (8) is the tensor product of crystals (the readers may regard this 
as just alternative notation). We call such elements of tensor products paths. 

The main ingredient of the formalism is the map called the combinatorial i?-matrices 

R : B^'' ® B^''' — > B^''' ® B^'' , . 

(a,6)®(c,d) ^ ic',d')(S){a',b'). 

In the present case A^^\ the explicit form of the map is 

a' = a + min(6, c) — min(a, d) 
b' = b — min(6, c) + min(a, d) 
c' = c — min(f), c) + min(a, d) 

d' = d + mm{b, c) — min(a, d). (3) 

An important point of the map R is that it has a deep mathematical origin as the inter- 
twining map that interchanges left and right of the tensor products of crystals. For the 
later purposes we introduce a vertex diagram for the map i? : a ® 6 H- 6' ® a' as follows: 



"'^In general, we can identify the Kirillov-Reshetikhin crystals B^''^ for type ^^^2^^ with the set of 
r X s semistandard tableaux with letters l,2,...,n. In this identification our (a, 6) is the height one 
semistandard tableau with a I's and b 2's. i?*"'^ corresponds to the Kirillov-Reshetikhin module naturally 
corresponding to the weight sA^ where is the r-th fundamental weight. 
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By a repeated use of the map R we define the time evolution of the box-ball systems 
Ti {I e Z>i) as follows. Let ui := (/, 0) be the empty box of capacity / and let b = 
bi ® b2 ® ■ ■ ■ ® bi be a given state of the box-ball system. We call ui the carrier. If 
necessary we put enough many empty boxes (1,0) on the right. Then we define 6']^, . . ., 
b'^ by the following diagram. 

bi 62 bt 



ui- 



■u 



(1) 



■u 
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Here the precise meaning of the diagram is as follows. We compute R : ui®bi ^ b'^® u] 
Then by using we compute R : uf^ ® 62 63 ®uf'\ We do this procedure recursively 
until the end of the state. Then we define 



Ti{b) ■.= b[®b'. 



(5) 



We can see that the time evolution rule given in Section [2] coincides with Too here. 

As a benefit of the definition by the crystal bases, we can show the quantum inte- 
grability of the box-ball system as the consequence of the Yang-Baxter relation for the 
combinatorial i?-matrices jFOYj . More precisely, we have 



Tin{b)=TkTi{b) 



(6) 



for arbitrary l,k G Z>i and states b. Moreover, we can construct conserved quantities of 
the box-ball system as follows. Let us define (see(jl])) 



(i-i) 



Eo{b) := 



(7) 
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ui and the energy function H : B^'^ ® B^''^' — > Z is defined by 

H{{a, b) ® (c, d)) := min(a, d). (8) 



Again an important point of the energy function is that it has a deep mathematical 
origin and is the consequence of the infinite dimensional symmetry of the quantum affine 
algebras. Let us consider the affinization of the crystal B 



Aff (5) = {b[d] \beB,deZ}. 



(9) 
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For elements of tensor products of Aff(i?), we introduce the afRne combinatorial R- 
matrices by 

Raff : bi[di] ® 62 M2] ^ - H{bi ® 62)] ® b[[di + H{bi ® ^2)], (10) 

where we have R : bi®b2 ^ 62®^! under the combinatorial R matrix. Then by the Yang- 
Baxter relation for the affine combinatorial i?-matrices we see that Ei are the conserved 
quantities of the box-ball systems [FOYj : 

Ei{n{b)) = Ei{b). (11) 

4 The rigged configurations 

Another important aspect of the box-ball systems is a connection with the rigged config- 
urations. In this section we give the definition of a special case of the rigged configuration 
bijection corresponding to the vector representation of type A^^\ Although this case is 
simpler than the general case, it is still nontrivial and we can see basic ideas of the theory. 

Originally the rigged configurations are discovered through insightful analysis of the 
Bethe ansatz for quantum integrable systems [KKR^ IKRj . The main ingredient of the 
theory is a bijection between the set of rigged configurations and elements of the tensor 
products of crystals. Such a bijection is generalized for highest weight elements of tensor 
products of the arbitrary Kirillov-Reshetikhin crystals of type An'^ and its mathematical 
theory is established by an important paper of Kirillov, Schilling and Shimozono jKSS] . 

In our case, a rigged configuration is composed of a Young diagram (called the con- 
figuration) and integers (called the riggings) associated with each row of the Young 
diagram. Let z/j (z = 1, . . . , (7) be the lengths of the rows of the configuration and let J, 
be the rigging associated with the row Vi. Then we represent the rigged configuration as 
(z^, J) = {{Pi, Ji)Yi=i- We call each {PiiJi) string. Although there is a characterization 
of the possible rigged configurations, we regard the set of the rigged configurations as the 
set of objects obtained by the map (in fact, bijection) 

$ : 6 I — >{v,J) (12) 

from arbitrary paths h. We call the bijection $ the rigged configuration bijection. 

Let us define the algorithm of the bijection $. For the given Young diagram v let 
Q^{v) be the number of boxes contained in the left I columns of p. Suppose that we are 
given the path b = &i ® 62 ® ■ ■ ■ 6l G (5^'^)®-^ where the positions of the balls bk = (0, 1) 
are given by = ki,k2,... from left to right. Let Pi{k^v) be the vacancy number 
defined by 

P,{k,v):=k-2Q,{v). (13) 

For example, we have ^2(16, yiH ) = 16 — 2 ■ 5 = 6. Suppose that a length L path b 
corresponds to the rigged configuration (z^, J). Then we call the string (z/j, Jj) singular if 
the rigging J, coincides with the corresponding vacancy number, that is, P,y^{L,iy) = Ji. 

The bijection $ is defined by a recursive procedure corresponding to the positions of 
balls /ci, ^2, . . .. We start from the empty rigged configuration. 
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1. Suppose that we have done the procedure up to kj^i and obtained the intermediate 
rigged configuration (77, /). 

2. For the next position kj, we do the following. Suppose that the rigged configuration 
(?7, /) corresponds to a length kj — 1 path. Compute the vacancy numbers Pr^X^j ~ 
1, rf) for all rows of r] and determine all the singular strings. 

3. If there is no singular string, add a length one row to the bottom of rj. Otherwise 
choose one of the longest singular string and add a box to the corresponding row. 
Denote by rj' the new configuration thus obtained|§ 

4. Define the new rigging /' as follows. For the strings that are not changed under 
f] —> f]', we choose the same riggings as before. Let rj'^ be the changed row under 
f] — > 1]'. Then define the new rigging by = Prj'.i^kj,!]') so that the string {r]'^,I-) is 
singular in [r]', I'). The output (?]', /') is the new rigged configuration corresponding 
to the length kj path. 

5. Repeat the same procedure for all kj. Let {u, J) be the final output. Then define 
$W = (^, J). 

A Mathematica package for the above procedure is available at |S3] . If we reverse all the 
procedure we obtain the algorithm for <l>~^. As examples, let us look at the example of 
the time evolution of the box-ball system at Section [2J In the first line, the positions of 
balls kj are 1,2,3,8. Then the computation of $ proceeds as follows: 



□ -1 
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(14) 



Here we put riggings on the right of the corresponding row and put kj above the corre- 
sponding arrows. Similarly, for the third line of the same example, we have 
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(15) 



5 The inverse scattering formalism 

The main observation on the relationship between the rigged configurations and the box- 
ball systems is that the rigged configuration bijection gives the inverse scattering formal- 
ism for the box-ball systems. In order to get the ideas of the result, let us compare the 
two examples in ( fT4l) and ( fT5l) . Then we see that the shapes of the configurations are 
same and the differences of the riggings are two times the lengths of the corresponding 
rows. Here we have the factor 2 in the change of riggings since we apply Too twice. 

In general, let b be the given state and let $(6) = {(z/j, Jj)}f^j^. Then we have |KOSTYj 

$(TKfe)) = {{i^i, J. + min(/, (16) 
^The order of rows is not essential in the definition of $. 
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This property is vahd for general box-ball systems including all cases that appeared in 
[HHIKTTl IFOYj . The proof of this fact heavily relies on a deep theorem of Kirillov- 
Schilling-Shimozono |KSS] ,^ Indeed, if we compare ( !T^ and f lT5|l we can see that this 
property is already nontrivial. To summarize, configurations are the conserved quantities 
(action variables) and the riggings are the linearlization parameter (angle variables) of 
the box-ball systems. Since $ is bijective, the rigged configurations give the complete set 
of the action and angle variables of the box-ball systemsjj 

Once we know that the rigged configurations are the underlying mathematical struc- 
ture of the box-ball systems, we can prove several fundamental properties of the box-ball 
systems. For example, the box-ball systems considered in |HHIKTTl IFOYj are shown to 
be solitonic by introducing a method to explicitly extract solitons from paths as elements 
of the affinization of the crystals |S1] . The main point of the proof of the result is to 
introduce a structure of the affine combinatorial i?-matrices on the rigged configurations 
via careful combinatorial arguments. We remark that the proof of the solitonic properties 
of the box-ball systems corresponding to the vector representation of type An^ is proved 
in |TNSj by taking certain ultradiscrete limit of an ordinary soliton system and an elegant 
alternative proof of their result is given in |FOY] by using the crystal bases. 

Another important problem that is solved by the rigged configuration bijection is 
the initial value problem of the box-ball systems |KSYlj . The result includes all the 
extensions considered in |HHIKTT| IFOYj . We note that an equivalent result for the 
case of the vector representation of A^^^ is rederived in |MIT2j . Let us explain the result 
for the case of the vector representation of A^i\ The main point is to give an explicit 
piecewise linear formula for the map : {ly, J) i — > b. For the given rigged configuration 
(z/, J) = {(z^j, Ji)}f^i, let us define the following ultradiscrete tau functions: 



Tr{k) := — min < '^^{Ji + — k)ni + min(^'j, v,j)ninj >, (r = 0, 1) 
" e{o,i}9 ^^^^^ J 



(17) 



where we denote n = (ni, . . . ,?T.g)EI Let us represent the k-th. element of the path b as 
bk = {1 — x{k), x{k)). Then we have the following analytic expression for the image b: 

x{k) = To{k) - Toik - 1) - ri(fc) + nik - 1). (18) 



■^If two tensor products b and b' are isomorphic under the combinatorial i?- matrices R : b ^ b' , we 
have 4>(6) = $(&') |KSS| Lemma 8.5]. The proof depends on a large part of the paper. 

''in fact if we restrict to consider the box-ball systems corresponding to the vector representation 
of type A^^\ we do not need to use heavy apparatus like rigged configurations. For example, jTTSj 
introduced a combinatorial method to obtain the conserved quantities. In jMITlj . a method to obtain 
the action and angle variables is derived, which is shown to be the special case of the rigged configurations 
jKSlj . In |KOTY| it is conjectured that the rigged configurations give the action and angle variables of 
the box-ball system corresponding to the vector representation of type A^^\ This problem is considered in 



|T1| with differently defined bijection. We remark that in [F] the Robinson-Schensted-Knuth algorithm 
is used to give some of the conserved quantities of the box-ball system corresponding to the vector 
representations of type An ^ (so called P-symbols arc conserved under the time evolutions). 

^If we consider paths with periodicities, these functions exactly coincide with the tropical Rie- 
mann theta function |KS2j . 
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Since the time evolution of the box-ball system is linearlized on the set of the rigged 
configurations, this result gives an explicit solution for the initial value problem of the 
box-ball systems. 

Sketch of the proof of fT^) . The main step of the proof is to show the following interpreta- 
tion of the tau functions. For the given path b = hi®h2®- ■ ■ , define Too(6) = b^i^^b^2^^' ' ' ; 
Tl{b) = 6^ ® fef ® ■ ■ ■ , and so on. Then we have to show the following interpretation: 

Tr{k) = (1 — r) X (number of balls in 61 ® 62 ® ■ ■ ■ ® fefc) 

+ ^(number of balls in b? ® O ■ ■ ■ (g) 6^*^). (19) 

i>l 

For example, in the example of Section [21 we have ro(8) = 9 and Ti(8) = 5. Since balls 
always move rightwards, the summation in the second term is always finite. From (fT9|) 
we can easily deduce ( !T8|) . 

Proof of (fT9!) proceeds as follows. From the expression (fTTj) we can construct determi- 
nants from which we obtain the tau functions as the ultradiscrete limit. Then by using 
a calculus of determinants we can show that the tau functions satisfy the ultradiscrete 
Hirota bilinear form. The Hirota bilinear form implies that the functions Tr corresponds 
to the same dynamics of the box-ball systems. Unfortunately this is not the whole story. 
The main difficulty is the fact that the analytic expression in ( ITTj) is very different from 
the combinatorial definition of the map and thus it is quite difficult to compare. 

We do this in the following way. The proof is induction on the rank n of An ■ Since 
we know that the tau functions satisfy the same dynamics of the box-ball systems, it is 
enough to consider a state T^{b) where N ^ 1. We call such a state the asymptotic 
state. Since we have the inverse scattering formalism which is the consequence of the 
most part of the paper |KSS] . we can easily obtain the corresponding asymptotic rigged 
configuration. Then we invoke the result of |S1] to reduce the problem to the case of A[j^2i 
(the case ^4^^'' can be shown by [SI]). This part is logically a bit complicated and we 
will omit the details. Here we remark that we use the fact that the tau functions for the 
general An^ have a similar recursive structure with respect to the rank and that we use 
the Yang-Baxter relations for the affine combinatorial /^-matrices to represent the right 
hand side of ( !T9|) by the energy function and the combinatorial /^-matrices. Thus the 
proof heavily utilizes the infinite dimensional symmetry behind the box-ball system. □ 

Finally we remark that the conserved quantities Ei of [FOYj indeed coincide with the 
rigged configurations [S2] : 

Ei{b) = Qi{u) (20) 

where $(6) = (z/, J)E There is a generalization of this formula for the most general rigged 

^In |T2| . Takagi introduced a scheme to factorize the dynamics of the box-bah systems of type A^^' 
into A^^' case by using the time evolution corresponding to the carrier of type B^'^. This scheme is 
rephrased into the rigged configuration language [KOSTY] Section 2.7] to factorize the map $ for general 
An"^ case into the map $ for A^^'^ case by using the B^'^ type time evolution. The proof of ([20]) in |S2| 
uses a refinement of the latter result. 
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configuration bijection of type An^ (see Section [7]). By fl2U]) we see tliat tlie configurations 
have a simple representation theoretical origin in terms of the combinatorial i?-matrices 
and the energy function. In fact, in [S2j it is shown that the quantities Ei (and suitable 
refinements) give enough information to reconstruct riggings. In this sense, the combina- 
torial algorithm of $ itself has a representation theoretical origin via the time evolutions 
of the box-ball systems. However the mathematical origin of the riggings is still unclear. 
In fact, the riggings depend on the information about the final positions where the cor- 
responding string is finally changed during the procedure $. This information is rather 
combinatorial and we cannot get rid of its difficulty even if we use the representation 
theoretical interpretation of the algorithm of $ discussed above. 

6 Interlude: the box-basket-ball systems 

In this section, we explain the basic properties of the box-basket-ball systems (BBBS 
for short) introduced by [LPS] . The starting point of the construction is to replace the 
combinatorial i?-matrices in the definition of the box-ball system by the whurl relations 
of |LP] . Rather non-trivially, the resulting dynamical system becomes a soliton system. 
The characteristic property of the BBBS is that the system contains the fermions (balls) 
and bosons (baskets) with mutual interaction between them. We remark here that the 
BBBS is different from the super-symmetric box-ball system of |HIj constructed from the 
crystals for the quantum superalgebra [BKKj since their system is the extension of the 
box-ball system by adding another kind of the fermionic particles. 

In order to obtain intuition about the model, it is convenient to start from a com- 
binatorial description of the time evolution of the BBBS. Let h = hi ® h2 ® ■ ■ ■ be the 
state of the BBBS. In this situation, each state is parametrized as a three dimensional 
vector (a, 6, c) & 1? . Our interpretation of each parameter is as follows; h is the number 
of baskets, c is the number of balls and a is the number of empty places that can fit 
extra balls. The meaning of such an interpretation will become clear when we explain the 
combinatorial description of the time evolution. 

In the rest of this section, we consider the following situation. We put many capacity 
one boxes on a line. If necessary, we put enough many empty boxes on the right of the 
state. As the rule, each box or basket can accommodate at most one ball whereas we 
can put more than one baskets on a box. Thus the balls are fermionic particles and the 
baskets are bosonic particles. There is a nontrivial interaction between the two kinds of 
particles by placing a ball within a basket. If necessary we assume that a ball is always 
placed in a box before placed in a basket. We introduce several definitions that will be 
used later. Let V = (1,0,0), F = (0,0,1), = (z + l,z,0) and Ui = where i > 1. 

Here we give several diagrams that represent these symbols: 

^ = U, ^ = IQI, B, = U, ^2=5. 

Now we explain the time evolution rule. We start from an initial state that contains 
finitely many baskets and balls. 
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Figure 1: Example of the time evolution of the BBBS. 

The time evolution of the BBBS: First, move every empty basket to the 
right one step. Full baskets are not moved. Second, consider each ball from 
left to right and move the ball to the next available empty box or basket. Each 
ball is moved exactly once. 

Note that if there is no basket, the above rule coincides with the one for the box-ball 
system. We give a simple but nontrivial example in Figure [1] 
The BBBS can be constructed from the whurl relation 

R : (a, b, c) ® {d, e, f) ^ {d' , e' , /') ® (a', b' , c') (21) 

where the explicit relations are 

a' = a — min(a + b,a + c,b + f) + min(e + c, d + c, d + b) 

b' = b — min(a + b,a + c,b + f) + min(a + e, d + f, e + f) 

c' = c — min(e + c, d + c, d + b) + min(a + e, d + f, e + f) 

d' = d + min(a + b,a + c,b + f) — min(e + c, d + c, d + b) 

e' = e + min(a + b,a + c,b + f) — min(a + e,d+ f,e + f) 

f' = f + min(e + c,d + c,d + b) - min(a + e,d + f,e + f). (22) 

We apply the whurl relation with ui := (/, 0, 0) to the diagram (jl]) to define the operators 
Ti {I G Z;>i). Then we can show that the above mentioned combinatorial definition of 
the time evolution coincides with T^o. Since the whurl relations satisfy the Yang-Baxter 
relation, we can show that TiT^^b) = T^Ti^b) for arbitrary /, /c e Z>i and states b. Thus 
the BBBS possesses the quantum integrability. We remark that since we do not know 
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the underlying symmetry of the whurl relations, we have not been able to construct a 
conserved quantity analogous to Ei. 

Even if we know that the BBBS is a quantum integrable system, it is far from clear 
whether the system is a soliton system. Below we explain that the system is indeed 
solitonic. For this purpose we classify solitary waves which do not change their shapes 
during the free propagations under Too- As the result, we see that there are the following 
two cases. 

1. A consecutive sequence of k balls Fk := FF ■ ■ ■ F. Under the free propagation by 
Too, Fk moves at velocity k. 

2. Any sequence of F,B,U which does not contain the consecutive subsequence FF 
or FU, which we call a slow soliton. Under the free propagation by Too, the slow 
solitons move at velocity 1. 

Note that F^ are the usual solitons of the box-ball system whereas the slow solitons are 
the new feature of the BBBS. 

Let us clarify what are the slow solitons. The answer comes from the analysis of the 
phase shift. Here the meaning of the phase shift is as follows. Let A and B are solitons 
on a line and suppose that they make collision during the time evolution and retain their 
original form after the collision. Then we compare the position of the soliton after the 
collision with the position of the corresponding soliton supposing that there is no collision. 
This difference (rightwards shift is positive) gives the phase shift. We summarize the basic 
physical properties of the fermionic solitons Fk and the bosonic solitons Ba^^Ba^ ■ ■ ■ Ba,^ in 
the following table. 





Fk 


Bai Ba2 ■ ■ ■ Ba^ 


velocity 


k 


1 


phase shift 


-2k 


-1 



Here the phase shift is defined by the scattering with Fi (/ > A;). For example, in the 
example in Section |2l we see that the length one soliton Fi get shifted by —2 after the 
collision with F^^ 

Let us look at two solitons Ti and Bi of velocity one. If we consider the scattering 
with Fk {k > 1), they get shifted by —2 and —1, respectively. To summarize, Ti and Bi 
have the same velocity whereas they have different values of the phase shift. Thus during 
the time evolutions we may have superposition of such states and this is the origin of the 
slow solitons. Therefore, in order to analyze the slow solitons, we make scatterings with 

^Let us mention the generalizations to the box-ball systems of the vector representations for types An'' . 
Then it is known that the phase shift coincides with the energy function (with a different normalization) 
between two solitons jFOYj . Here we identify freely propagating solitons with the semistandard tableaux 
and regard them as the elements of crystals B^^'^ of types Note that since we are neglecting 

all I's (empty places), we have A^^}_-^ here. In |Slj . it is generalized to include all cases considered in 
[HHIKTTl IFOY] and the scatterings of soitons are identified with the affine combinatorial i?-matrices 
([T0|) where each soliton corresponds to the truncated rigged configurations. 
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many F^s and decompose them into elementary solitons Fi and Bi. For example, the 
example in Figure [1] shows the decomposition of the slow soliton U2 into two elementary 
solitons Fi and -82- 

Based on these observations, we define the solitons of the states of the BBBS as the 
elementary solitons Fi and which we can obtain by scattering with many additional 
Ffc's. Let us define the amplitudes of Fi and Bi by / and i, respectively. Then we can show 
that the number and amplitudes of the solitons are preserved during the time evolution of 
the BBBS. Moreover we can show that scatterings of multiple solitons can be decomposed 
into two body scatterings. Hence we see that the BBBS is solitonic. 

7 Generalizations and further developments 

In the most of the present note, we only think about the simplest possible case, namely 
the vector representation of type A^-^^ We do so in order to provide the basic ideas 
without getting into the technical complexities. In fact, one of the nice features of our 
approach is its universality. For example, the definition of the box-ball system in (jl]) has 
straightforward generalizations for the tensor products of the Kirillov-Reshetikhin crystals 
B^'^ for the quantum affine algebras of types other than An^ . Here B'^''^ is the Kirillov- 
Reshetikhin crystals corresponding to the weight sA^, where is the r-th fundamental 
weight. In this case, instead of using ui in we use the classically highest weight 
element of B^''^. Then we denote by T^'^ the resulting time evolutions. Again the box and 
ball interpretation of the time evolution provides a nice way to get intuition about the 
generalized models. For example, for the box-ball systems corresponding to the vector 
representations of general non-exceptional affine algebras, there is an interpretation of 
2^1,00 terms of particles and anti-particles with pair creations/annihilations jHKT] . In 
this final section, we will give comments on the methods of the generalizations and further 
properties. 

Known extensions. The rigged configuration is known to have many extensions. In- 
deed it is expected that such a bijection exists for the arbitrary Kirillov-Reshetikhin 
crystals corresponding to general affine quantum algebras. As mentioned in Section HI 
the bijection for type An^ is already constructed in full generalities. Apart from this case, 
we have the following generalizations. 

• ^ B^'^ for arbitrary non-exceptional affine algebras |0SS2] . 

• (8)^5^-1 for type D^n^ [84]. 

• (g).^!'^' for type d'^^ 

• B''' of type D'n^ jOSST] . 

• (g)Si'i for type E^^^ [0^2] . 
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We remark that the combinatorial algorithms involved in these extensions share many 
common features and the philosophy which underlies these extensions is the same. We 
also remark that all these results are related with the highest weight elements of tensor 
products of crystals. However, if we think about the box-ball systems we encounter the 
rigged configurations for not necessarily highest weight elements. This extension is quite 
natural. Indeed the algorithm presented in Section |4] does apply to both cases without any 
change. Therefore it is quite natural to consider the Kashiwara operators (analogue of the 
Chevalley generators in the crystals setting) on the set of the rigged configurations. This 
is achieved in |S5] for all simply laced cases. Remarkably, the definition of the Kashiwara 
operators for the all cases considered in [S5] is uniform. 

The method of generalizations. As examples of the generalizations, let us consider 
the cases An or Dn ^ . Then the rigged configurations take the following form: 

J) = ((/^W, J^'^), {^^'\ J^'^), ■ ■ ■ , (^("\ /"^)) (23) 

together with the Young diagrams which is determined by the shape of the tensor 
product B = i?''"*' by the following rule: each B^''^^ in B corresponds to the length Si 
row of Note that we should consider that each (//'•"^ J'-"^) corresponds to the node 
a G Jo of the Dynkin diagram (without the 0-node, see |Klj ) for the affine algebras. For 
example, the vacancy number for the present case takes the following form 

Pt\i^) = Qiifi^"^) - 2Q,{u^'^^) + Q^(^^'^)^ (24) 

bSilo, b^a 

where a ^ b means that the nodes a and h are connected by a single edge on the Dynkin 
diagram. Note that the definition ( !T3|) corresponds to the special case a = 1 and (jS^^ = 
{1^). There is a nice characterization of the highest weight rigged configurations. For the 
given rigged configuration {ly, J), if 

P?,(^)> Jf^>0 (25) 

is satisfied by all the strings Ji""^), then (z/, J) is the highest weight rigged configu- 
ration. 

The algorithm $ is almost parallel to the definition given in Section H] (see, for exam- 
ple, [S2l Appendix A] for details). To get a feeling of the algorithm, suppose that we have 
a letter a in a type B^'^ path (in Section HI we described the case a = 2). Here we iden- 
tify the elements of crystals B^''^ with semistandard tableaux or Kashiwara-Nakashima 
tableaux (generalizations of the semistandard tableau, see [KNj ). Then we have to add 
a box to each of u^°'~^\ iy^'^~'^\ ■ ■ ■ , u^^^ in this order. The rule for the addition to u^"-'^^ 
is exactly the same one given in Section |H Suppose that we have added a box to the 
column of Then we look for the longest singular string of (z/'^'^~^\ J*^""^-*) 

whose length does not exceed to determine where to add a box. We do this recur- 

sively until (i^*-^-*, J^^^) by recursively defining £*^^^'s. Finally change the riggings by using 
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the new vacancy numbers as in Section HI For the modifications required for the negative 
letters in type Dn \ see, for example, [OSSlj . Roughly speaking, we do an almost similar 
procedure twice (for a, first proceed from u^""^ to u^^^ and next to the left from u^"'^ as 
above) following the crystal graph for the vector representations B^'^ of type Dn\ 

The inverse scattering formalism (fT6ll holds almost identically for the general cases. 
For arbitrary T''' of type A^n^ |K()STY] and for T^'' of type dI^^ |KSY2] . it is known 
that the only change caused by T'"'* is the shift in the rigging 

minis, u^^)) (26) 

and all the other places do not change. We expect that a parallel formalism should exist 
for all types of the quantum affine algebras once the corresponding rigged configuration 
bijection is established. 

Similarly, the relation (!20|) has the following straightforward generalization for ar- 
bitrary rigged configurations of type An^ (not necessarily highest weight). In (jl]) and 
([7]), we use the time evolution T*"'* instead of T;(= T^''). Then we can define (6) as 
generalizations of Ei{b). Then we have |S2] 

E'''ib) = Qsiiy^'''^). (27) 

On the other hand, the ultradiscrete tau functions formalism is only available for the 
case 5^'^* of type An^ . Perhaps we need to thoroughly understand the dynamics of 
the box-ball systems for general cases (say, the case B^^'^ of type An^). 

Further properties. So far we have explained that the rigged configurations behave 
very nicely with respect to the box-ball systems. In particular, the rigged configurations 
have the concrete mathematical meaning as the action and angle variables for the box-ball 
systems. As the final remarks we explain there are equally remarkable properties of the 
rigged configurations with respect to other mathematical problems. In most cases, the 
rigged configurations behave surprisingly simply with respect to global and deep structures 
of the corresponding algebras which are usually difficult to realize. 

• The combinatorial i?-matrices become trivial on the level of the rigged configura- 
tions. If the two tensor products are isomorphic under the combinatorial i?-matrices 
i? : 6 H- 6', we have $(6) = $(fe'). Remind that the combinatorial i?-matrices for 
general situations are highly complicated objects. This property is confirmed in all 
known cases and we expect that it is true for arbitrary quantum affine algebras. 

• The Schiitzenberger involution and its generalizations become almost trivial oper- 
ation (see, for example, |KSSt ISS]): we take complements of all the riggings with 
respect to the corresponding vacancy numbers. Note that in this case we consider 
only the highest weight rigged configurations which satisfy (123]) . 

• In |0S1] a new kind of bijection \Ef for the rigged configurations is introduced. The 
map \E' gives one to one correspondence for the following two sets; (i) the set of the 
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highest weight rigged configurations for arbitrary non-exceptional quantum affine 
algebras of sufficiently large rank, and (ii) the set of pairs of the highest weight 
rigged configurations of type An^ and the Littlewood-Richardson tableaux. In ex- 
periments, we can see that the map \1/ coincides with (and generalizes) the global 
involution exchanging the nodes and n of the Dynkin diagram of the algebra (if 
such an involution exists, see jLOS] ). Remarkably, the construction of the algo- 
rithm is quite simple and does not depend on the choices of the corresponding 
non-exceptional algebras. Indeed, the algorithm coincides with the type An^ rigged 
configuration bijection if we change left and right in the definition (as if we are using 
mirrors). The Littlewood-Richardson tableaux naturally appear as the recording 
tableaux. We remark that such a correspondence is very difficult to construct if we 
do not use the rigged configurations (see [S7]). 

• In |S6] the affine Kashiwara operators for type Dn"^ are realized via the Dynkin 
involution exchanging the nodes and 1. The realization relies on a rather non- 
trivial bijection between the Kashiwara-Nakashima tableaux and a combinatorial 
objects called the plus-minus diagrams. Then the involution is realized as changing 
columns of the plus-minus diagrams. In jOSSlj we see that the plus-minus diagrams 
essentially coincide with the rigged configurations. Thus we can realize the Dynkin 
involution 1 as a transformation on the rigged configurations. 

However the main point of the result is not the practical values. Rather, the result 
reveals that the crystal structure of the corresponding case is essentially governed 
by the rigged configurations. Note that in this case the rigged configurations for 
non-highest weight elements play the role. 

Concluding words. We have seen that the rigged configurations have very special 
properties which are usually difficult to see so that it is tempting to say that they are one 
of the canonical realizations of the Kirillov-Reshetikhin crystals. Not only they give a 
nice presentation, they also have concrete mathematical meanings and it seems that they 
originate from deep aspects of the infinite dimensional symmetry of the quantum affine 
algebras. Although the theory of the rigged configurations is still in a very early stage, 
we expect that the progress of the theory will give unique insights into the nature of the 
symmetry of the quantum affine algebras. 
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